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ABSTRACT 

Eclipsing binaries (EBs) provide critical laboratories for empirically testing 
predictions of theoretical models of stellar structure and evolution. Pre-main 
sequence (PMS) EBs are particularly valuable, both due to their rarity and the 
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highly dynamic nature of PMS evolution, such that a dense grid of PMS EBs is 
required to properly calibrate theoretical PMS models. Analyzing multi-epoch, 
multi-color light curves for ~2400 candidate Orion Nebula Cluster (ONC) mem- 
bers from our Warm Spitzer Exploration Science Program YSOVAR, we have 
identified twelve stars whose light curves show eclipse features. Four of these 12 
EBs are previously known. Supplementing our light curves with follow-up opti- 
cal and near-infrared spectroscopy, we establish two of the candidates as likely 
field EBs lying behind the ONC. We confirm the remaining six candidate sys- 
tems, however, as newly identified ONC PMS EBs. These systems increase the 
number of known PMS EBs by over 50%, and include the highest mass {6^ Ori 
E, for which we provide a complete set of well determined parameters including 
component masses of 2.807 and 2.797 Mq) and longest period (ISOY J053505.71- 
052354.1, P~20 days) PMS EBs currently known. In two cases {6^ Ori E and 
ISOY J053526. 88-044730. 7), enough photometric and spectroscopic data exists to 
attempt an orbit solution and derive the system parameters. For the remaining 
systems, we combine our data with literature information to provide a prelim- 
inary characterization sufficient to guide follow-up investigations of these rare, 
benchmark systems. 

Subject headings: open clusters and associations: individual (Orion) — stars: pre- 
main sequence — stars: binaries: eclipsing — stars: variables: general 



1. Introduction 



The Orion Nebula Cluster (ONC) contains several thousand members, and since it is 
nearby, it provides an excellent empirical laboratory to study the physical properties of pre- 
main sequence (PMS) stars and brown dwarfs. The ONC is particularly useful for comparison 
of the obser ved luminosities and effective temper atures of PMS stars to theoretical model 
predictions ( iHillenbrandl 119971 : iDa Rio et al.ll2010[ ). Such comparisons can, in theory, allow 
an estimate of the masses of individual stars as well as both the mean age and the age spread 
for the stars in a cluster. For such estimates to be meaningful, however, the theoretical tracks 
and isochrones must be vetted against observations to insure that they are well calibrated. 
Empirical measurements of the masses, radii, and temperatures of stars, over a range of 
masses, are necessary for understanding stellar evolution and for deriving well-calibrated 
theoretical models. 



The most rigorous means to measure precise stellar properties is via identification and 
characterization of PMS eclipsing binary (EB) stars because, through a complete analysis of 
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spectroscopy and photometry of these systems, the individual masses, radii, temperatures, 
and absolute luminosities of the two stars can be accurately derived. However, the identifica- 
tion of such systems is difficult due to the need for monitoring observations and the fact that 
the system must have an inclination close to 90° to be detected. The paucity of known PMS 
EEs has meant that the theoretical models lack rigorous empirical confirmation, and thus 
that masses derived from those tracks have significant systematic uncertainties associated 
with them ( iHillenbrand fc Whitd l2004j ) . The situation has begun to change recently with 
the advent of sensitive wide-field cameras, robotic telescopes and automated photometry 
pipelines, allowing deep, wide, long duration time series monitoring programs to be con- 
ducted. These programs have now led to the identification o f seven low-mass PM S EEs with 
i ndividual masses lo wer than 1.5 M r;,: RXJ 0529.4+0041A JCovino et al.lboool). VI 174 O ri 



fIStassun et al.ll2004f). 2MJ05 3 5-05 fIStassun et alJ l2006l. 120071 ). JW 3 80 fllrwin et al 
Par 1802 JCardle et al.lboosi Istassun et al.lboosi i. ASAS J0528-F03 dstempels et al 



2007) 



20081 ) 



and MML 53 dlfebb et al.lboioi I2OI1I ). These systems have components ranging in mass from 
0.036 Mo (2MJ0535-05E) to 1.38 Mq (ASAS J0528+03A). All but ASAS J0528+03, MML 
53, and RXJ 0529.4-1-0041 are associated with the ONC. There are just another four EE sys- 
tems where either onl y the secondar y is on the PMS or both comp onents are more massive 



PMS stars: EK Cep f jPopperl Il987h . TY CrA fICasev et al.lll998h . RS Cha f lAlecian et al. 



2OO51 . 120071 ) and perhaps also V578 Mon ( IGarcia et al.ll201ll ) although its components are 
E type stars that may already be on the main sequence. T he masses of a handful PM S 
low-mass sta r s hav e been measured usi n g other methods (see IHillenbrand fc Whitd ( 120041 ) ; 
Eoden et al.l (|2005[ ) ; ISimon et al.l ( 120001 ) ; iTognelli et al.l ( 120111 ) for a summary on dynamical 
mass determination and calibration of PMS tracks). 

Despite the recent discoveries, it is still important to search for additional PMS EEs. 
Main sequence solar-type stars are well described by state-of-the-art stellar evolution mod- 
els (i.e., observations agree well with theoretical isochrones); however, recent measurements 
of the stellar properties of low-mass dwarfs and young PMS stars remain problematic for 
the existing models. Recent work has suggested that in addition to mass and age, other 
parameters, such as magnetic field strength or rotation, may be necessary to fully charac- 
terize young, low mass stars. Magnetic fields in young, rapidly rotating low mass stars are 
thought to inhibit convection and thereby cause those stars to have l arger radii and cooler 
temperatures than would otherwise be the case ( jMorales et al.ll2010t iMacdonald fc MuUan 
20101). This effect has been invoked to explain the properties of the 1 Myr old brown dwarf 
EE in Orion (2 M0535-05) where the more massive cornponent is unexpectedly cooler than 



its co mpanion ( iReiners et al.ll2007t iChabrier et al.l 120071 : iStassun et al.ll2007t iMohanty et al. 



2OIOI ). The effects of magnetic fiel ds on stellar structu re are not included in the models 
and are not completely understood ( IChabrier et al.ll2007l ). Most PMS eclipsing binaries dis- 
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covered to date have short enough orbital periods that their components are expected to 
have their rotation periods tidally locked to the orbital period. Therefore, they are likely 
to be rapidl y rota ting and have strong magnetic fields, and hence have infiated radii (see 
Kraus et al.l ( 1201 ll )). Identification of PMS EBs with longer periods, where tidal locking is 
not expected, could offer a direct test of the proposed link between magnetic fields, rotation 
and radii. Moreover, 2M0535-05 remains the only known brown dwarf EB. 

In this paper, we report the identification and initial characterization of six new PMS 
EB candidates in the ONC, discovered as part of the YSQVAR (Young Stell ar Object VARi- 
ability) Spitzer Exploration Science program ( iMorales-Calderon et al.ll201l[ ). These systems 
have been overlooked in the past probably because the area surrounding the Trapezium 
stars is filled with bright nebulosity, making the optical photometry very unreliable, but 
also because cadence and duration of observations in some previous studies were not ideal. 
Details of the discovery observations are reported in Section 2. Follow-up observations are 
presented in Section 3 and 4. In section 5 we provide a description of the available data and 
preliminary analysis for the 6 PMS EBs, based in most cases upon the light curves, in order 
to provide a basic initial characterization. Two additional systems are likely field star EBs 
lying behind the ONC and are described in the Appendix. 



2. Discovery Light Curves: The 2009 Campaign 
2.1. Spitzer Observations in 2009 



As part of the Spitzer Explor ation Science progra m YSOVAR, we used about 250 hours 
of warm Spitzer Space Telescope (IWerner et al.l 12004) observing time to monitor ~0.9 deg^ 



of the ONC with the Infrared Array Camera (IRAC; iFazio et al.ll2004l ) at 3.6 and 4.5 fim in 
Fall 2009. The observed area was broken into five segments (five astronomical observation 
requests or AORs) with a central region of ~20'x25' and four fianking fields. The central 
part was observed in full array mode with 1.2 seconds exposure time at 20 dithering positions 
to avoid saturation by the bright nebulosity around the Trapezium stars. The remaining four 
segments of the map were observed in High Dynamic Range mode with exposure times of 10.4 
and 0.4 seconds, at 4 dithering positions. These observations were taken f or 40 days in Fall 
2009, with ~2 epochs each day. We used the IDL package Cluster Grinder ( ICutermuth et al. 



20091 ) which, starting from the basic calibrated data (BCD) released by the Spitzer Science 
Center, builds the combined mosaic for each AOR at each epoch and performs aperture 
photometry on the mosaics. Because we dithered while mapping and then extracted the 
photometry from the mosaiced image, each star at each epoch will have contributions to its 
image from four to as many as 80 BCDs; the time and magnitudes we report for each epoch 
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correspond to the average for all of those frames. 



From these data, we constructed light curves for ~10000 point sources in our FOV. 
About 2400 of these stars were identified as probable ONC members, including ~1400 stars 
with probable mid-IR excesses (Megeath et al. 2012) and ~1000 additional stars identi- 



and variability studies ( 


^arenaeo 


1954; 


Jones & Walker 


1988; 


McNamara 


1976 


1996; 


Tobin et al. 


2009; 


Carpenter et al. 


2001 




Getman et al. 


2005|), but lackin 



significant enough to be included in the Megeath et al. catalog of YSO candidates. Most 
of these latter stars should be weak-lined T Tauris (WTTs), though some may be sources 
with excesses that escaped previous detection. While our primary goal was to use these 
data to investigate the structure of the inner disk and time-variable accretion in YSOs 
with circumstellar disks, these observations also provided a treasure of data on all types 
of PMS variability. Fu r ther d iscussion of this dataset and initial results can be found in 
Morales- Calderon et all J201l[ l. 



In order to search for eclipsing systems, we first ran a Box-Least-Squares (BLS; lKovacs et al 



20021 ) periodogram on all 2400 of these sources and then automatically scanned the folded 
light curves looking for the signature of detached EEs - short duration flux dips that are 
present and similar in both IRAC channels and light curves that are otherwise approximately 
constant. In addition, we also visually inspected the light curves searching for signatures 
that our code could have missed. We also examined in a similar fashion the light curves of 
the ~1000 anonymous stars in our catalog brighter than [4.5] = 13.5 mag (which corresponds 
to the quiescent magnitude of our faintest EB candidate). 

Fro m this process, we i dentified 12 EB candidates. Fou r of them, 2MASS 05352184' 



Stassun et al. 



0546085 fIStassun et al.ll2006h . V11 74 Ori fIStassun et al.ll2004h . Parl802 fICargile et al.ll2008 



20081 ). and JW 380 (llrwin et al.ll2007n are already known PMS EBs. Another 



six systems ar e newly identified ONC PMS EBs -includin g 6^ Ori E, a known PMS spectro- 
scopic binary (jCostero et al.ll2006l ; iHerbig fc GriffinI 120061 ) which was fiagged as potentially 
eclipsing, but no confirmatory photometry had been obtained until now- and are shown in 
detail in the following sections. The remaining two sources are not known Orion members 
from the literature and, based on our analysis of the available data, we believe are likely to 
be background stars, though still echpsing binaries (see Appendix Al and A2). All but one 
(ISOY J053454.31-045413.0, hereafter: ISOY J0534-0454) of the new PMS EB seem to have 
no detectable IR excess and they are likely to be WTTs. 

We used implementation of the BLS period-finding algorithm available at the NExScI 



site (http://exoplanetarchive.ipac.caltech.edu) to determine if the fiux dips in each light 
curve are indeed periodic and to identify the best period. In some cases, when we had 
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several bandpasses or epochs well-separated in time we tweaked the obtained period by 
hand to obtain a more accurate one consistent with all epochs and bandpasses. For the 
four previously known systems, our derived periods agree, within errors, with the published 
period. Our new PMS EB candidates are listed in Table [H The first column of Table 
[1] displays the name of the source which is formed by an acronym (ISOY: Initial Spitzer 
Orion YSO) followed by the coordinates of the source. For the sake of simplicity we will 
use reduced names (i.e.: ISOY JHHMM-DDMM) throughout the text and figures except for 
ISOY J053515. 76-052309. 9, for which we will use its most common name in the literature: 
9"^ Ori E. The remaining columns in Tabled] show available broad band photometry, spectral 
types, and periods. Finding charts from the Spitzer 4.5 /im images for our 6 PMS EBs are 
provided in Figure [1] 



2.2. Ground-based Observations in 2009 



To complement our Spitzer data, we obtained contemporaneous Ic and J photometry, 
usually for smaller areas within the Spitzer mosaic. The main source of J-band monitoring 
data was the United Kingdom Infrared Telescope Wide Field Camera (UKIRT/WFCAM). 
For the Ic band, the New Mexico State University/ Apache Point Observatory (NMSU/APO) 
Im telescope and NasaCam at the 31" telescope at Lowell Observatory contributed the bulk 
of the monitoring data. We performed differential aperture photometry on the ground-based 
data. In each dataset, for each object, an artificial reference level was created by selecting 
30 nearby isolated stars; we iteratively eliminated those with larger photometric scatter or 
evidence of variability. Zero-points for the photometry were established by reference to data 
for non-variable stars in our FOV from the Two Micron All-Sky Survey (2MASS) P oint 
Source Catalog (jCutri et al.ll2003l ) for J, and by reference to data in iHillenbrandl (119971 ) for 



All of the 2009 time series data are publicly available through the YSOVAR database 



(http://ysovar.ipac.caltech.edu/). The 2009 light curves for the six new PMS EB are shown 
in Figure [2l Where we have photometry at the right epoch, the ground-based data corrob- 
orate the Spitzer eclipses. 



3. Follow-up Light Curves: The 2010-2011 Campaign 



To confirm the identification of the new EBs and to refine measurements of their orbital 
periods and eclipse shapes, we obtained additional photometry in 2010 and 2011. These 



- 8 - 



observations are described in the next sections and summarized in Table |2j 



3.1. Spitzer Observations in 2010 and 2011 

Four of our new PMS EB candidates: 9^ Ori E, ISOY J0535-0525, ISOY J0535-0523, 
and ISOY J0535-0522, are located within 3' of 6^ Ori C, the brightest star in the Trapezium 
cluster. These stars were observed in IRAC full-frame mapping mode with a frametime of 
2 seconds and 20 small dithers, at a cadence of 7 to 8 times per day, for the period Nov. 
19-29, 2010. Another two systems, ISOY J0535-0447 and ISOY J0536-0500, were observed 
for just two days in November (Nov. 20-22) 2010 at the same cadence, but with a 12 second 
frametime and four small dithers. The observations were designed to detect at least one 
eclipse (either primary or secondary) in each system. In addition, sparser data (~8 epochs 
in 30 days) was also obtained for all of our EB candidates in Fall 2010 with the same setup 
as in the 2009 mapping observations (same spatial coverage and integration times). Finally, 
in November 2011 we obtained three observations per day for ISOY J0536-0500 and ISOY 
J0534-0454 for about 20 days with a 12 second frametime and four dithers at each epoch. 
The IRAC light curves for 2010 and 2011 were extracted in the same way as for the 2009 
data. Figure |3] shows these light curves for the 6 PMS EBs. 

For ISOY J0535-0447, we also obtained ~10 hrs of IRAC staring mode observations 
at [4.5] on 2011 April 27, scheduled to cover the secondary eclipse of ISOY J0535-0447. 
In this case we simply stared at the target object, with no dithering, in order to provide 
more accurate relative photometry. A total of 2713 consecutive frames of data with 10.4 
sec integration time were taken. The integration time was selected so that the number 
of electrons would be high; however, to make sure that the target would not saturate, we 
placed it in the center of four pixels. For the light curve extraction of these data, we 
performed aperture photometry on the corrected BCDs (cBCDs) supplied by the Spitzer 
Science Center. The cBCDs are calibrated frames with empirical corrections for artifacts 
due to b right sources. However, instrumental effects such as the pixel phase effect can still 



remain (IMorales-Calderon et al.l l2006l : ICody fc Hillenbrandl l201l[ ). We inspected our light 



curve searching for signs of this instrumental signature (flux strongly correlated with sub- 
pixel position) but, if present, its effect on the light curve is negligible compared to the 
amplitude of the observed variations. 
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3.2. Ground-Based Observations in 2010 and 2011 

As in 2009, we also obtained contemporaneous Jc and J pliotometry in 2010 and 2011. 
The data were obtained from the UKIRT/WFCAM, the NMSU/APO Im telescope, the 
Wide Field Imager on the 2.2 m telescope at La Silla Observatory, and the 40" telescope at 
the United States Naval Observatory (USNO) Flagstaff Station. The J band UKIRT mosaic 
covers the whole area observed by Spitzer, and we obtained one observation per night for 
about 30 days in Fall 2010 (between Oct. 29 and Dec. 15). For 10 of those nights two 
observations per night were obtained. APO, WFI and USNO data were obtained for the 
Trapezium area with a cadence of once a night between October 7 and December 14 for 
APO, between November 21 and November 28 for WFI, and for several hours a night for 
34 nights between October 30 and December 28 for USNO. Finally, the largest effort was 
performed on the Northern area where ISOY J0535-0447 and ISOY J0534-0449 are located. 
We used the USNO 40" telescope to monitor a ~23'x23' region for several hours per night 
for 38 nights in 2010 and 35 nights in 2011 (a total of 15137 frames were obtained). The 
light curves were extracted in a similar fashion as for the 2009 data. 

The light curves for our 6 PMS EBs with the follow-up photometry from 2010 and 2011 
can be seen in Figures |3] and HI All the YSOVAR time series photometry at all available 
bandpasses for the 6 PMS EBs can be found in Tables E] to [TTJ 



4. Spectroscopic Follow-up 
4.1. Low- Resolution Spectroscopy 



We obtained moderate-resolution near-infrared spectra of three of our PMS EE cand i- 
dates using Triplespec at the Palomar 5m telescope (IWilson et al.l 120041 iHerter et al.ll2008l ). 
TripleSpec covers the wavelength range from 1 to 2.45 /im simultaneously at a spectral res- 
olution of ~2700. The entrance slit is 1x30 arcsecond and the spectrum is spread over 
five orders. Observations were obtained at two slit positions, using a standard ABBA nod 
sequence. Spectra of ISOY J0535-0447, ISOY J0535-0522, and ISOY J0535-0523 with total 
exposure times of 720, 720 and 1200 seconds, respectively, were obtained on September 20 
and 21, 2011; another 720 second spectrum of ISOY J0535-0447 was obtained on November 
27, 2010 with the same setup. 

We obtained similar spectra for ISOY J0535-0454 and I SOY J0536-0500 us ing the Triple- 
Spec spectrograph mounted on the APO ARC3.5 telescope (IWilson et al.ll2004j ). APO Triple- 
Spec instrument is a twin of Palomar TripleSpec and thus covers the same wavelength range 
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(1 to 2.45 /iin) at similar resolution. Both sources were observed on October 30, 2011 at two 
slit positions and with a total exposure time of 5280 seconds. The 1x45 arcsecond slit was 
used yielding a resolution of ~2300. 

The spectra were reduced using a custom version of the IDL-based SpeXTool reduction 
package (jCushing et al.ll2004j ). as modified to process data obtained by the Palomar and 
APO TripleSpec instruments. Sky subtraction was performed by differencing spectral frames 
obtained with the target at the A and B slit positions; visual inspection of the sky subtracted 
frames identified apparent sky line artifacts in the sky subtracted 2-dimensional spectra of 
ISOY J0535-0523 (primarily the He I line at 1.083 microns) and numerous lines throughout 
the JHK bands for ISOY J0535-0522. These residual 'sky lines' are due to line emission 
from the substantial nebulous regions near J0535-0523 and J0535-0522, which can be easily 
seen in the 4.5 /xm finding charts presented in Figure [TJ 

Each target's flattened, sky-subtracted, extracted, and wavelength calibrated spectrum 
was corrected for telluric absor ption and flux cali brated using the Interactive Data Language 
(IDL) XTELLCOR package (jVacca et al.l |200J) and observations of AOV stars obtained 
nearby in time {5t < 30 min) and airmass {5 sec z < 0.2) to each target observation. 

We used these spectra to derive spe ctral types for our targets by comparison with 
dwarf sta ndards in the in the S pex library (IRayner et al.ll2009l ) and young TWA and Taurus 
members (ILuhman et al.ll2006[ ). The derived spectral types are presented in Table [Hand the 
K band region of the spectra are shown in Figure O 



4.2. High-Resolution Spectroscopy 

On December 13, 2010 we obtained single epoch high resolution spectra for ISOY J0535- 
0447, ISOY J 0535-0522, and IS OY J0535-0523, using the High Resolution Echelle Spectrom- 



eter (HIRES; IVogt et al.lll994l ) at the Keck I telescope. In all cases we used the red cross- 
disperser and a 0.8" width slit providing a resolution of R=A/AA =50,000 and covering the 
wavelength range from 4310 to 8360 A. The exposure times were 120, 300, and 600 seconds 
respectively. In addition, for ISOY J0535-0447, we took an extra epoch on March 15, 2011 
and on December 10, 2011 we took single epoch spectra for ISOY J0535-0525 and ISOY 
J0534-0449 with exposure times of 900 and 1800 respectively, keeping the same setup. 

Data reduction was performed with "MAKEE" (MAuna Kea Echelle ExtractiorB), the 
standard data reduction package for the HIRES instrument which, starting from the raw im- 



* http: / /www2 .keck.hawaii.edu/inst /hires /makeewww / 



- 11 - 



ages, produces wavelength calibrated spectra. A section of the spectra, showing the lithium 
absorption, for the 4 PMS EBs observed with HIRES can be seen in Figure [61 

For ISOY J0535-0522, ISOY J0535-0525, and I SOY J0534-0454 we also ob tained spectra 



with the NIRSPEC infrared echelle spectrograph (IMcLean et al.lll998l . l2000f ) at Keck II on 
November 9, 2011, January 11, 2012, and January 15, 2012 respectively. We used the H band, 
centered at 1.555/im, with the 0.288x24" slit, yielding a resolution of ~30,000. The total 
integration time was 8-16 minutes. Data re duction was accomplished with the REDSPEC 
packageQ. Additional details are provided in IpratJ hmi^ . 



For ISOY J0535-0447, we have tried to obtain enough data to attempt to derive an 
orbit solution. For this purpose we obtained one additional high resolution spectrum using 
the echelle spectrograph at the Mayall 4m telescope at Kitt Peak National Observatory 
(KPNO) on September 24, 2010, and five epochs between November 22, 2010 and February 
5, 2011 with NIRSPEC/Keck II . Mayall 4m observations were obtained using the Echelle 
spectrograph with the 'long red' camera, to cover 4540 - 7620 A at an average resolution of 
32,000. We obtained two exposures of 900 seconds each. A standard echelle-data reduction 
was applied. This includes background subtraction, cosmic-hit removal, flat-fielding and 
wavelength calibration. NIRSPEC observations were obtained using the setup described 
above and total integration times were 6 to 8 minutes. 

We measured heliocentric radial velocities by cross correlation of the target spectrum 
with those of standard stars of similar spectral type using several orders. The accuracy of 
the radial velocities is 0.9-2 km/s, as measured by spectra of bright F and G stars of known 
radial velocity obtained on the same nights as the targets' spectra. The derived radial velocity 
measurements from each spectroscopic observation are listed in Table |3l ISOY J0535-0525 
and ISOY J0534-0454 showed clear double lines in their spectra. Both the H band and the 
optical spectrum of ISOY J0535-0522 are suggestive of a blend of two spectra but additional 
observations will be required to confirm the presence of double lines. None of the spectra of 
ISOY J0535-0447 showed double lines (see Sec. 5.2 for further discussion). 



5. Results and Discussion 

To put or EE systems into context we show in Fig [7] a J vs. J- [3.6] color-magnitude 
diagram showing the location of the new PMS EBs (blue circles; Ori E is not shown due 
to the lack of shorter wavelength data) and the previously known ones (red squares) together 



t 



http://www2.keck.hawaii.edu/inst/nirspec/redspec/index.html 
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with the Orion known members. 

In this section we combine Spitzer hght curve data with follow-up ground-based pho- 
tometry and radial velocity measurements, as well as data from the literature, to perform 
an initial analysis of the physical parameters for each of the six PMS EBs discovered in 
the YSOVAR program. In cases where we have the most data, including well-sampled light 
curves at multiple wavelengths and radial velocity measurements, we attempt to measure 
the component stellar masses, radii, temperatures, and other parameters. For most sources 
we have only the YSOVAR light curves, and data based on colors and/or spectral types, 
so our derived properties are fairly rudimentary. We present these preliminary analyses to 
guide follow-up investigations of these rare, benchmark systems. For the remaining two EB 
systems, our analysis of the available data suggests that while eclipsing, the systems are not 
young members of Orion. We summarize the properties of these two systems in Appendix 
A. 



5.1. Ori E 



6^ Ori E, the " fifth" member of the Orion Trapeziurn, is a known double-lined spectro- 
scopic binary (SB2; Iflerbig fc Gr"iffi^l2nn6l : ICostero et ailbood . boosi ). ICostero et alT ( |2008l ) 
found that the two stars have nearly identical spectra, with spectral types of GO IV to 
G5 III (with a most likely spectral type of G2 IV). Strong Li I A6708 absorption and mod- 
erate Ca II K emission are observed, confirming the pre-main-sequence status of the binary 
and, hence, its membership in the Orion Nebula Cluster. Given the likely large stellar radii 
and relat ively tight orbit of the binary, it has been suspected that 6^ Ori E might undergo 
eclipses (IHerbig fc Griffin! |2006| ). However, its proximity to the brighter Trapezium star 
6^ Ori A has in the past made it difficult to obtain reliable photometry. The high spatial 
resolution and improved brightness contrast in the infrared afforded by Spitzer makes an 
investigation of the EB nature of 9^ Ori E possible for the first time. 

We show in Figures [2] and [3] (middle right panels) the Spitzer light curves obtained by 
us for 6^ Ori E in 2009 and 2010, respectively. Both light curves show periodic eclipses in 
both the [3.6] and [4.5] IRAC bands, showing for the first time that 6^ Ori E is indeed an 
EB. The eclipses are quite shallow, indicating that this is a grazing eclipse. 

In addition to the eclipses, there are indications of variations in the overall fiux of the 
system. The median level of the Spitzer light curves is ~64 mmag fainter in 2009 compared 
to 2010, whereas the [4.5] fiux observed by us in 2010 is the same as was observed by Megeath 
et al. (2012) in 2004. We have explored the possibility that this small change in brightness 
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is due to the photometry being affected by the bright nebulosity surrounding the Trapezium 
stars by checking the light curves of nearby sources. None of them showed a similar change in 
flux from one season to the next and thus we believe that the brightness change observed is 
real. Indeed, 6^ Ori E has also be en reported to ex hibit intrinsic variation in previous studies 
at X-ray and radio wavelengths. IKu et al.l ( 1l982l ) noted that the X-ray flux fluctuated on a 
timescale of a few ks, and the nearly continuous X-ray light curve spanning 13 d obtained 
in 2003 by the Chandra Orion Ultradeep Project (COUP; IStelzer et al.l l2005l ) showed a 
gradual 20% brightening of the star followed by an abrupt spike at about twice the original 
brightness. The s tar was also detected in Very Large Array observations at 2 and 6 cm by 
(IFelli et al.lll993al Jbl). during which variations of ~50% in radio flux were seen. These X-ray 
and radio variations may not be surprising given the SB2 nature of the system, however the 
cause of the infrared variations remains unclear. The most likely cause of these intrinsic flux 
variations is stellar activity (star spots; flares; other coronal structures). Here we simply 
note their existence and adjust the 2009 IRAC light curve to match the continuum flux of 
the 2010 data in order to permit an initial modeling of the infrared light curves. The time 
series data are provided in Table [61 



We combine these light curve data with the radial velocitie s measured by 



Costero et al. 



( 2008 ) to produce a complete analysis of the binary. We used the Wilson fc DevinnevI package 



(1971 , updated 2005; hereafter WD) as implemented in the PHOEBE code by lPrsa &: Z witter 
( 120051 ). The code uses Kurucz stellar atmosphere models to represent the unde rlying stellar 
f luxes across the bandpasses, as well as the empirical limb darkening laws of Ivan Hamme 
( 119931 ). The radial velocities and light curve data are fit simultaneously and self-consistently. 
The WD code does not include filter profiles or limb darkening coefficients for the IRAC band- 
passes, so for this initial modeling, we approximated the [3.6] bandpass with the Johnson 
L bandpass and the [4.5] bandpass with the Johnson M bandpass. Strictly speaking, these 
differences between the bandpasses used for the observations and for the modeling introduce 
a systematic error in the model fluxes, but these differences should be small given the ap- 
proximate similarity and broadband nature of the bandpasses, as well as the survey-grade 
quality of the light curves used here. The orbital period and time of periastron passag e for 
the system were held flxed at the values previously determined by ICostero et al.l ( 120081 ) . 



The IRAC light curve flts resulting from our WD modeling are shown in Figure [SI and 
the derived system parameters summarized in Table [H The precision and limited wavelength 
coverage of our light curve data do not permit highly precise determination of the stellar 
radii or temperatures, however we are able to provide constraints on the sum of the stellar 
radii and on the ratio of the components' effective temperatures. These parameters and 
the system inclination show some degeneracy (Figure [9l), which we incorporate into our 
parameter uncertainties in Table [H The mass ratio of the system is extremely well determined 
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from the excellent quality radial velocity data of ICostero et al.l (120081 ) . for which we find 
q = 0.9965 ± 0.0065; the two stars are of very nearly identical mass. The system inclination 
is now also very well constrained from our light curve fits (i~74°), and thus we are able to 
report individual masses for 9^ Ori E of Mi=2.807 M© and M2=2.797 Mq, accurate to ~2%. 
It is thus the highest mass EB yet discovered with clear PMS nature. 

The temperature ratio that we determine for 6^ Ori E is T2/Ti=1.12 ± 0.08. This is 
roughly consistent with the expectation of identical temperatures if the stars are of truly 
identical mass; however, the measured mass ratio also allows for a small mass difference which 
would then also accommodate a small temperature difference. Higher precision light curves 
over a larger range of wavelengths will be required to more firmly pin down the temperature 
ratio. The inferred temperature ratio is also modestly dependent upon the adopted absolute 
temperatures of the stars. Here we have assumed Ti = 6000 K based on the reported spectral 
type, but this too needs to be improved. 

The sum of the radii that we measure is Ri + R2 = 12.5 ± 0.6 Kq. It is premature 
to attempt a detailed comparison to the predictions of pre-main-sequence stellar evolution 
models, however these r adii appear t o be b roadly consistent with expectations. For example, 
the model isochrones of ISiess et al.l ( 120001 ) predict a radius sum of 11.1-14.5 R© for ages in 
the range 0.3-2 Myr. 



5.2. ISOY J0535-0447 



ISOY J0535-0447 is a known proper- motion member of the ONC (ITian et al.l Il996l ). 
and our own spectroscopy shows strong lithium absorption (see Figure E]), further confirming 
ISOY J0535-0447 as a PMS member of the ONC. From our TripleSpec data (see Figure [5]) 
we derive a s pectral type of K0-K 2 and extinction of Ay=1.5-2 based on a comparison to 
SpeX dwarfs (IRayner et al.lll989l ). Our HIRES spectrum yields a spectral type of Kl (± 
1 subtype) based on comparison with dwarf standards. This star has also been included 
in previous variability surveys of the ONC, but has not he r etofor e been identified as an 
EB. Neither the near -IR variability survey of ICarpenter et al.l (120011 ) nor the Chandra X-ray 
variability survey o f iRamirez et al.l (|2004j) found ISOY J0535-0447 to be photometrically 



variable. Similarly, iTobin et al. 
for SBs in the ONC. 



(I2OO9I ) did not flag this star in their radial-velocity survey 



That the EB nature of this star was previously missed both photometrically and spec- 
troscopically is likely due to a confluence of several factors: the orbital period is very close 
to an integer number of days (P ~ 3.9 d; see below), the optical echpses are relatively shal- 
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low (see Figure 131 lower right panel), the secondary eclipse is extremely shallow and only 
readily detectable in the IRAC bands, and the optical spectrum does not reveal a secondary 
spectrum but instead appears as a single star of spectral type ~K1. 

The secondary eclipses in this system are sufficiently shallow that from our discovery 
Spitzer light curves it was not immediately evident whether secondary eclipses were present 
at all. In addition, our intensive follow-up light curves from the ground in Jc-band, which 
are of excellent precision, exhibit out-of-eclipse variations whose amplitude is larger than the 
secondary eclipse and which vary with a period slightly different from the orbital period (see 
below), making ready identification of the secondary eclipse difficult. Therefore we obtained 
a follow-up IRAC light curve at the expected secondary eclipse time (see Figure S]) definitively 
proving the existence of a faint secondary star, and providing a very strong constraint for 
the modeling of the system, which we now describe. 

The most comprehensive light curve data were obtained with the USNO 40" telescope at 
Ic band, with high cadence sampling over several hours per night on about 70 nights in 2010 
and 2011. Given the strong out-of-eclipse variations present in the /c-band data which might 
dominate over the eclipse signal with Fourier based period-search methods, we determined 
the eclipse period using a BLS period-finding algorithm and then we manually adjusted the 
resultant best period by coni b ining the data from the different bandpasses and adding the 



monitoring data from iRebullI ( 120011 ). We determine the following system ephemeris, which 



we adopt throughout our analysis: 



HJDo=2455126.26 
P=3.905625±0.000030 days 



As mentioned, the 1^ band light curve exhibits a strong periodic modulation superposed 
on the eclipses. In addition, a flare is visible around HJD 2455502 (day 375 in Fig. |3]). 
These features strongly indicate the presence of spots on the primary star modulated on 
its rotation period (the secondary is very likely too faint to produce such a strong spot 
modulation signal). We find this spot modulation to have a period slightly longer than the 
EB period (4.001 d vs. 3.906 d), indicating a modest departure from full synchronization 
of the system, and which produces the stroboscopic effect that causes the apparent secular 
changes in the primary eclipse depths over time (Figure |3]). In addition, we found that the 
spot modulation is not well fit by a single sinusoid, but rather appears to have two sinusoidal 
components that moreover evolve modestly in their relative amplitudes and phases. This 
is presumably due to changes in the spot coverage of the stellar surface causing shape and 
amplitude changes in the out-of-eclipse modulations. 
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In principle it should be possible to model the system light curves incorporating a full 
physical spot model. For this initial investigation, we have performed a rectification of the 
/c-band light curve in an attempt to remove the complex out-of-eclipse variations. We fit the 
light curve with a two-component Fourier model, with a single principal period of 4.001 d 
(see above). To allow for the possibility that the variability changed in amplitude and/or 
phase on long timescales, we fit and subtracted the model to the 2010 and 2011 /c-band 
data separately. Overall, the light curve fits are reasonably good (see Figure [10]). Some 
unmodelled variation in the J^-band light curve still remains and the residuals affect the 
primary eclipse slightly and perhaps also the secondary eclipse depth. The model for the 
primary eclipse at longer wavelengths in Figure [10] appears by eye to overestimate the eclipse 
depth. Due to the uncertainties and limited sampling in the photometry for the primary 
eclipse at these wavelengths, relative to the higher-precision, higher- cadence staring mode 
data of the secondary eclipse, our fitting routine places more weight on the accuracy of the 
fit to the secondary eclipse. The primary eclipse is well fit in our shortest wavelength, and 
thus it is unlikely that there are large systematic errors in the secondary and/or primary 
radii. However, there may be some small errors in the assumed extinction-corrected colors 
of the components in this system. 

None of the spectra taken for this EB candidate show double lines; however, they do 
show small radial velocity (RV) variations (see Table [3]). The small RV amplitude together 
with the very shallow secondary eclipse depth point to the companion being much lower 
mass than the primary. 

We have used the few available radial velocity measurements, together with the light 
curve data, for a first attempt to obtain the physical parameters of the system. The pho- 
tometric times of minima indicated no sign of orbital eccentricity, which is expected for a 
short-period binary such as this, so we assumed circular orbits in the light curve analysis. 
In addition, since we only have RV measurements for the primary, we had to make one more 
assumption for the orbital solution. Thus, we adopted a total mass of ~0.9 M© (by fixing the 
semi-major axis, a=lOR0) and then fit for the mass ratio and the center-of-mass velocity. In 
Table [5] we present the orbital parameters of ISOY J0535-0447 resulting from our best orbit 
solution and display this solution along with the radial- velocity measurements in Figures [TU] 
and[TT] Note that our derived center-of-mass velocity for ISOY J0535-0447 of 7 = 30.4 km/s 
i s slightly larger than the currently accepted ONC systemic cluster velocity of 25± 2 km/s 



( jSicilia-Aguilar et al.ll2005l ). The presence of spots on the stella r surfaces can also int roduce 



distortions that affect the observed radial velocities; however, IStassun et al.l (l2004j ) found 



that the spot effect on the RVs is ~1 km/s and probably negligible for our data. 



With estimated masses for the system components of 0.83 Mq and 0.05 M©, the sec- 
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ondary seems to be in the brown dwarf regime, filling an important gap in the current 
census of PMS EBs. Note however that this is a preliminary result; the inferred temperature 
of the secondary is slightly warmer than what would be expected for a brown dwarf. We 
have an ongoing study of this system to further characterize the physical properties of the 
components of ISOY J0535-0447 and measure their masses more accurately. 



5.3. ISOY J0535-0522 



ISOY J0535-0522 is located at about 1' northeast of the Trapezium. It has been included 
in several previous studies, but it is in a very complicated region with bright structured 
nebulosity which has probably prevented its identification as an EB before. ISOY J0535- 
0522 was detecte d as an X-ray sourc e and hence confirmed as a probable member of the 



ONC by CO UP flGetman et al. 



fiagged possible fiaring star by 
young stars in the Orion Nebula. 



20051) . It was also included in the near-IR variability study 



performed by lCarpenter et al.l (12001) where it was not fou nd to be variable, however, it was 



Feigelson et al.l (120021 ) in their study of X-ray-emitting 



ISOY J0535-0522 was classified as a late K or early M star by iLuhman et al.l ( 120001 ). 
however, our TSpec data (see i^-band in Figure^]) is very similar to that of ISOY J0535-0447, 
but with higher extinction (At,=11.5-12), suggesting a spectral type near KO. The hydrogen 
lines appear contaminated by nebular emission and without those lines a K3 spectral type 
cannot be ruled out. A spectral type of K0-K3 and large Ay are more consistent with the 
location of the system in a color-magnitude diagram (see Figure [7]) than the previously 
reported spectral type. A Keck/HIRES single epoch spectrum yields a spectral type of KO 
± 2 subtypes and shows Call core emission in the triplet lines and lithium absorption (see 
Figure |6]), both features confirming its PMS status. Both the NIRSPEC and HIRES spectra 
of ISOY J0535-0522 are suggestive of the SB2 nature of the system (see the HIRES spectrum 
in Figure [6]) and the IRAC [3.6] and [4.5], and Ic band light curves from 2010 confirm the 
eclipses discovered in 2009. 

The light curves of this EB candidate are presented in Figures [2] and |3] (lower left 
panels) for the 2009 and 2010 observing runs respectively and the time series are available 
in Table [HI The 2009 light curves show a trend in the [4.5] data between HJD 2455147 and 
2455152 (between days 20 and 25 in Fig. 12]) and between HJD 2455157 and 2455161 (between 
days 30 and 34 in Fig. |3]) that is not followed by the [3.6] data and we believe is not real. 
Those datapoints were excluded from the following analysis. The 2010 data did not show 
similar artifacts but did show an ascending trend that lasted for all the 2010 period. Again, 
the [3.6] data did not show the same behavior and, given that the nebulosity is brightest at 
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[4.5], we believe that our photometry was affected by it. We corrected the 2010 IRAC [4.5] 
time series by subtracting a hnear fit to the whole raw light curve in 2010 (the light curve 
shown in Figure [3] as well as the time series provided in Table [H] are already corrected from 
that effect). 

The light curves include over 10 clearly discernible eclipses which have allowed us to 
perform a periodogram analysis using the BLS period-finding algorithm. The resulting 
ephemeris is: 



HJDo=2455128.79 
P=5.6175± 0.0015 days 



Using these parameters, we present the folded light curves in Figure [121 A rough 
comparison of the primary and secondary eclipse depth yields a temperature ratio of T2/T1 ~ 
0.94 which, assuming a temperature of 5250 K for the KO primary (based on the PMS 
temperature scale of iKenyon &: Hartmanrull995l ). produces a temperature of ~4935 K for 
the secondary, which would then be a K2 star. 



5.4. ISOY J0535-0523 



ISOY J0535-0523 was labeled as a proper motion member of the ONC by I Jones &: Walker 



(119881 ) and it is located close to the Trapezium s tars as well (<3^ to the wes t of the Trapez 



i um). It has bee n labeled as non- varia ble by both [Carpenter et al.l (j200l[ ) and lFeigelson et al. 



(120021 ). In 1994, IStassun et al.l (119991 ) performed an Ic band monitoring of the ONC; a sim 
ilar study was performed by the Monitor Project during 2004 to 2006 (jlrwin et al.l 120071 ) 
neither group identified the star as being an EB. 



A spectral type of M4.5 has been previously reported (lHillenbrandlll997l ) and our TSpec 
data (see Figure [5]) show that it is indeed a mid-M star, M5-M5.5, with A^=1.5 based on 
a comparison of steam bands to members of TWA and Taurus. This is in agreement with 
the location of ISOY J0535-0523 in the color-magnitude diagram showed in Figure [71 Our 
single epoch HIRES spectrum yields an spectral type of M6 ± 1 subtype and shows lithium 
absorption as can be seen in Figure [61 The spectrum does not show obvious double lines. 

The light curves for 2009 and 2010 are shown in Figu res [21 and [31 where a handful 
of eclipses are noticeable. We have added the dataset from [Stassun et al.[ (119991 ) and the 
Monitor Project to our 1^ band observations in order to use the long time baseline to improve 
the accuracy of our period analysis. From these light curves we have derived a tentative 
ephemeris as: 
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HJDo=2455136.7 
Period=20.485± O.Sdays. 

However, among all the data, only one eclipse was well mapped. This and the fact that the 
eccentricity is not exactly zero makes it more difficult to pin down the period. Thus, we 
provide above the most likely period and note that the error in the period determination 
for this source is larger than for the other systems. The folded light curve is provided in 
Figure [T3l The Monitor Project Jc band data (from 2004) con&m the eclipse s ignatures seen 



i n our data, however there is some disagreement with the /c-band data from IStassun et al. 



( I1999I ) (at phase~-0.4) probably due to a small error in the period accumulated over the 
years. Due to the scattering of the photometry in the eclipses, it is not straight forward 
to determine which one is the primary or the secondary and thus an estimation of the 
temperature ratio is not possible. 

ISOY J0535-0523 has a period far longer than any of the known PMS EEs and much 
longer than the expected rotational period. It may thus be the only known EB that is 
probably not spun up by tidal effects. If dynamo activity increases the radii of most EB 
stars, ISOY 0535-0523 could be the only system known whose components' radii are not 
inflated by these tidal effects, thus allowing its properties to be better compared to the 
current generation of theoretical models providing a direct test of the link between magnetic 
fields, rotation, and radii. Performing this test will require an RV solution in order to fully 
characterize the system. 



5.5. ISOY J0534-0454 



ISOY J0534-0454 is catalogued as having a disk based on its mid-IR excess (Megeath et 
al. 2012). Its spectral energy distribution (SED) can be seen in Figure UM ISOY J0534-0454 
is the only EB candi date that we found arn ong the disked obje c ts. It has been labeled as 
non variable both by lCarpenter et al.l (l200l[ ) and lRamirez et al.l ( 120041 ). 



ISOY J0534-0454 is heavily extincted and so we lack optical photometry. Our TSpec 
data indicates a mid-M spectral type, which is consistent with its location on a J vs. J-[3.6] 
color- magnitude diagram (see Figure [7]) assuming high extinction. The Keck/NIRSPEC data 
shows double lines in the spectrum indicating that the system is an SB2. The light curve of 
ISOY J0534-0454 in 2009 (see Figure [2]) shows a handful of eclipse features clearly seen at 
least at IRAC wavelengths (which has a higher cadence than the J UKIRT data). For this 
system, we do not have high cadence observations in 2010 (Fig. [3]) and the lower cadence 
data (~12 datapoints in 40 days) does not show any indication of the presence of an eclipse. 
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However the Spitzer data from 2011, where a few echpses are present, further confirm the 
EB nature of the system and yield the ephemeris as: 

HJDo=2555127.56 
Period=5.1993± 0.0004 days. 

The folded light curve using that period can be seen in Figure [151 We note that the light 
curve is remarkably clean despite the presence of a disk. ISOY J0534-0454 seems to be a 
system with e ~0 and with two similar components with T2/T1 ~ 0.96 bas ed on the eclipse 



depth. Assuming a temperature of 3125 K for an M5 primary using the iLuhmanI (119991 ) 
temperature scale for young M stars yields a temperature of ~2990 K for the secondary, 
and thus an M6 star. The secondary may then be a brown dwarf given its spectral type (the 
known brown dwarf EB in Orion is an M6.5). 



5.6. ISOY J0535-0525 



ISOY J0535-0525 was selected as a proper motion member of the ONC by 



( 1988 ). It has been labeled a long term variable by iFeigelson et al.l (120021 ): iHerbst et al. 



Jones fc Walker 



( I2OOOI ) found that it was a periodic variable with a period of ~5.7 days, which is consistent 
with our photometry. A wide r ange of spectral ty pes are assigned to this target in the liter- 
ature: MO. 4, K4, and K7-M1 ( IHillenbrandl 1 19971 ) . probably due to its location close to the 
Trapezium cluster and the variable extincti on in the re g ion. I t was included on the radial- 
velocity survey for SBs in the ONC done by lTobin et al.l ( 120091 ) where it was not flagged as a 
probable binary (Vrad=22.9±2.25). However, our high resolution spectra of this source shows 
the presence of double lines and thus confirms its binary nature. Furthermore, the HIRES 
spectrum we obtained shows it to have a strong lithium absorption feature, confirming its 
youth (see Figure |6]). Thus, we can classify it as a PMS EB. 

However, only our 2009 IRAC data show evidence for this designation (see Figures [2] 
and [3]). The 2009 IRAC light curve shows a well-defined approximately sinusoidal variation 
presumably reflective of rotational modulation of a spotted star. However, at two epochs - 
HJD 2455132 and 2455156 (days 5 and 29 in FigM) - the photometry of both IRAC channels 
show a 0.15 mag dip of duration less than one day. We do not have shorter- wavelength data 
at the times of these dips. If the dips are due to eclipses, and if no eclipses were missed in 
the IRAC data, then the period would be of order 24 days. There is another event at HJD 
2455528 (day 401 in Fig. [3]) and several other single point deviants that we cannot confirm 
as eclipse features. 
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This could be another very useful PMS EB if its period were indeed confirmed to be 
~24 days, since it would then be the longest period PMS EB known and could serve as an 
additional empirical test of PMS isochrones for objects that are not tidally locked to short 
rotation periods. 

6. Conclusions 

We have identified six new EB star systems that are beheved to be members of the 
Orion Nebula Cluster and thus arc PMS EBs. For one of them, 9^ Ori E, we provide an 
orbital solution with a complete set of well determined parameters. The masses derived 
for its components are 2.807 and 2.797 Mq and thus 6^ Ori E is the most massive EB 
known to date with clear PMS nature. For a second system, ISOY J0535-0447, we provide 
a preliminary orbital solution based on our light curves and available radial velocity data. 
For the other four, we provide periods and relatively well-defined light curves, often based 
on multi-year data. One of these systems, ISOY J0535-0523, is the longest period (P~20 
days) PMS EB currently known. Such a long period suggests that the components will be 
less affected by tidal effects providing a direct test between magnetic fields, rotation and 
radii. Another two systems seem to have secondary components in the brown dwarf domain 
providing, if confirmed, additional examples to the only one known PMS brown dwarf EB. 

These systems increase the number of known PMS EBs by over 50%, and the unique 
properties of several of these systems ensure that they will offer the potential for considerably 
improving the empirical calibration of the PMS models for low mass stars and brown dwarfs. 
However, additional observations are needed to fully characterize the systems, in particular 
high resolution spectroscopy at several epochs is needed to derive the orbital parameters of 
these systems. 
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A. Additional Objects of Potential Interest 

We have identified two other objects which seem to be eclipsing binaries, but where the 
data we collected suggest that the objects are not bona-fide Orion members. We describe 
these two objects below. 



A.l. ISOY J0536-0500 



ISOY J0536-0500 has not been previously labeled as an ONC member. It is fainter than 
all previously known EBs and, because it is located further to the NE of the Trapezium, 
it has not been included in most of the studies of the region. Its brightness and colors are 
very similar to those of the known brown-dwarf EE 2M0535-05 as is shown in Figure [71 Its 
position in the color-magnitude diagram indicated that ISOY 0536-0500 could comprise two 
substellar objects of even lower mass than those in 2M0535-05. 

ISOY J0536-0500 light curves show several eclipses, both in IRAC [3.6] a nd [4.5] and J 



bands are clear. This source was also included in the Ic band monitoring of IStassun et al. 



(119991 ) and we have used their 1^ band data to obtain a more accurate period. The timing 
analysis gives a period of 3.5705 days and the ephemeris we calculate is: 



HJDo=2455128.74 
Period=3. 570535 ± 0.00007 days 



The phased light curve based on these parameters can be seen in Figure Again, the 
If. band data from 1994 confirms the eclipse features. 

The problem with this picture is that the near-IR spectrum which we obtained does not 
correspond to that of a late-type dwarf. Instead, the TSpec spectrum indicates the star is 
much hotter, with an inferred spectral type of G2 or earlier (see Figure E]). The broad-band 
colors are inconsistent with this spectral type unless the star has a reddening of order A^, ~ 
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6. The location of the star in a CMD (see Figure [7]) is inconsistent with this reddening unless 
the object is not a member of the ONC. Given the properties of this star, we suggest that 
it is a field EB located behind the Orion molecular cloud and seen through an edge of the 
cloud. 



A.2. ISOY J0534-0449 



ISOY J0534-0449 was identified by lStassun et al.l (119991 ) as the fast est rotator in a stud y 
of the period distribution in the ONC (P=0.27 days, Star 1161 of IStassun et al.lll999l ). 
Rebulll (l200l[ ) phot ometrically monitored t he outer ONC finding the same period (Star 1450 
in their study), and lCarpenter et al.l ( 1200 ll ) finds our EB candidate periodic as well but with 
a period of 3.14 days (their cadence did not allow them to find shorter periods). 

Figur e [T71 includes a ll the Y SOVAR data from 2009 to 2011 plu s the /^ - band monitoring 
data from IStassun et al.l (119991 ) (obtained on December 1994) and iRebullI (120011 ) (obtained 
between December 1995 and February 1997) folded with a period of 0.5424. The light curve 
shape appears remarkably similar in all these datasets, spanning 17 years at multiple epochs. 



Its pubhshed spectral type is K5 (lRebullll200ll )): however, its location in a J vs. J— [3.6] 
diagram shows it to be much fainter than the other Orion members of similar spectral type - 
but not significantly reddened (see Figure [7]). In addition, our single epoch HIRES spectrum 
shows that about the only feature visibile is a fairly strong Hq, absorption line instead of in 
emission as it would be expected for a young mid-K star. If it were a single, field rapid rotator, 
presumably the light variations would be due to large, non-axisymmetrically distributed 
spots. However, it seems very unlikely to us that such spots could yield such a stable light 
curve shape over such a long period. Instead, we believe the object is more likely to be a 
field contact binary composed of two ~G dwarfs. 
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Fig. 1. — Finding charts from the Spitzcr 4.5;Lim images for each of our 6 PMS EEs ordered 
by right ascension. The size of each image stamp is 3'x3'. Note that the intensity scale 
varies from one image to the next. 
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Fig. 2. — Light curves for the 6 new Orion PMS EBs discovered by the YSOVAR program 
from 2009 data. The symbols are •: [3.6]; o: [4.5]; *: J (UKIRT); and +: (APO & 
LOWELL). [3.6] and [4.5] magnitudes have been plotted in the right and left vertical axis 
respectively, /c and J light curves, if present, have been shifted in the y-axis to match the 
mean IRAC values. Mean colors are stated in each panel. Note that the x-axis is different 
depending on the available data. 
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Fig. 3. — Light curves for the 6 new Orion PMS EBs discovered by the YSOVAR program 
from the follow up observations in 2010 and 2011. The symbols are •: [3.6]; o: [4.5]; *: J 
(UKIRT); and +: Ic (APO). [3.6] and [4.5] magnitudes have been plotted in the right and 
left vertical axis respectively. Ic and J light curves, if present, have been shifted in the y-axis 
to match the mean IRAC values. Mean colors are stated in each panel. Note that the x-axis 
is different depending on the available data. 
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Fig. 4. — Follow-up IRAC [4.5] staring mode light curve for ISOY J0535-0447 obtained on 
April 27, 2011 showing the secondary echpse. 
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Fig. 5.— TripleSpec K band spectra for ISOY J0536-0500, ISOY J0535-0447, ISOY J0535- 
0522, ISOY J0534-0454, and ISOY J0535-0523 from top to bottom. The spectra have been 
normahzed and offset in the y axis for clarity. The derived spectral types for each source are 
shown. The Bracket gamma and 2.12 micron H2 lines in the spectrum of ISOY J0535-0522 
and ISOY J0535-0523 respectively, are likely due to nebular emission. 
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Fig. 6. — HIRES spectra for the four PMS EB for which we have high resolution optical 
spectra showing the presence of lithium absorption. From top to bottom: J0535-0523, J0535- 
0525, J0535-0522, and J0535-0447. The spectra have been normalized and offset in the y 
axis for clarity. 
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Fig. 7. — J vs. J-[3.6] color-magnitude diagram showing the location of our 6 PMS EBs (blue 
circles) and the previously known ones (red squares) together with known Orion members 
(grey circles). The position of the two suspected field EBs is also shown with green triangles. 
A 1 Myr isochrone is plotted as a solid line. 6^ Ori E is not shown due to the lack of J-band 
data. 





Fig. 8. — Light curves for 6^ Ori E phased with a period of 9.89520 days in both IRAC 
channels, 3.6 /im (top) and 4.5 fim (bottom). The best fit model light curve is produced 
(see Sec 5.1) using the parameters from Table [Hand it is overplotted as a solid red line. The 
residuals of the fit can be seen in the lower part of each panel. 
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Fig. 9. — Dependence of the radii (left) and the temperature ratio (right) on the inchnation 
for 6^ Ori E showing the degeneracy of the fitting process. The results of our fit are marked 
with a dot and the contours represent 1, 2, and 3 sigma confidence levels. Inclination and 
combined radius are highly degenerated; inclination and temperature ratio are less so, but 
also provide weaker constraints. 
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Fig. 10. — Light curves for ISOY J0535-0447 folded with a period of 3.905625 days. From 
top to bottom: Jc-band (rectified as described in Sec. 5.2), J-band, 3.6 fim, and 4.5 fim. The 
best fit model light curves are produced for each filter using the parameters from Table [5] 
and are overplotted as solid red lines. Note that the IRAC [4.5] staring data around the 
secondary eclipse provides a very strong constraint for the modeling of the system. The 
residuals of the fit can be seen in the lower part of each panel. Even around the secondary 
eclipse the fit is not perfect due to the stellar variability. 
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Fig. 11. — Radial velocity measurements for the primary star of ISOY J0535-0447. A model 
radial velocity curve is generated using the parameters from Table [5] and is overplotted with 
a solid red line. The residuals of the fit are shown in the bottom of the panel and have an 
rms of 1.07 km/s. 
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Fig. 12. — Phased light curve for ISOY J0535-0522 folded with a period of 5.6175 days where 
all the data from 2009 and 2010 are plotted. The symbols are •: [3.6] IRAC; o: [4.5] IRAC; 
and +: Ic WFI;. [3.6] and Ic light curves have been shifted in the y-axis to match the mean 
[4.5] values ([3.6]-[4.5]=0.31 mag, /c-[4.5]=6.65 mag). 
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Fig. 13. — Phased light curve for ISOY J0535-0523 folded with a period of 20.485 days 
where all the data from 2009 and 2010 are p l otted . The symbols are •: [3.6] IRAC; o: [4.5] 
IRAC; and +: APO+WFI+ lstassun et aD f ll999h + the Monitor Project. [3.6] and light 
curve have been shifted in the y-axis to match the mean [4.5] value ([3.6]-[4.5]=0.05, Jc- 
[4.5] =3. 15). Note that there are some Ic band data around phase=-0.4 that d oesn't line up. 
Our p hotometry phases well with the Monitor Project data but the data from lStassun et al. 
(1l999h shows some disagreement probably due to noisy photometry or a small error in the 
period accumulated over the years (see Sec. 5.4) 
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Fig. 14. — Spectral energy distribution for ISOY J0534-0454 showing the presence of a 
moderate mid-IR excess indicative of a circumstellar disk. The combined photometry for 
the system is used. The data come from 2MASS (diamonds) and IRAC (black circles for 
data coming from Megeath et al. 2012 and grey circles being the median values from our 
YSOVAR light curves). For comparison a reddened M5 Kurucz model normalized at K band 
is shown. 
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Fig. 15. — Phased light curves for ISOY J0534-0454 folded with a period of 5.1995 days 
where all the data from 2009 to 2011 are plotted. The symbols are •: [3.6] IRAC; o: [4.5] 
IRAC; and *: J UKIRT. [3.6] and J light curves have been shifted in the y-axis to match 
the mean [4.5] value ([3.6]-[4.5]=0.21, J-[4.5]=4.49). 
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Fig. 16. — Phased light curves for ISOY J0536-0500 folded with a period of 3.5705 days 
where all the data from 2009 and 2010 ar e plotted. Th e symb ols are •: [3.6] IRAC; o: 
[4.5] IRAC; *: J UKIRT, +: J, band from Istassun et all Jl999[ ). [3.6] , J, h light curves 
have been shifted in the y-axis to match the mean [4.5] value ([3.6]-[4.5]=0.02, J-[4.5]=1.44, 
I-[4.5]=3.81). 
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Fig. 17.— Phased light curves for ISOY J0534-0449 folded with a period of 0.542413 days 
where all our data fr om 2 009 and 2010 are plotted together with 1994-1997 Ic band data from 



Stassun et al. 


(1999 


) and 


RebuU 


(2001) 



. The symbols are •: [3.6] IRAC; o: [4.5] IRAC; *: J 



UKIRT; and +: 1^ USNO+SMMV+ROl. [3.6] , J, and h band light curves have been shifted 
in the y-axis to match the mean [4.5] value ([3.6]-[4.5]=0.02, J-[4.5]=1.41, /c-[4.5]=2.81). 



Table 1. New eclipsing binary candidates 



Source^ 


Most Common Namc*^ 


V (mag) 


Ic (mag) 


,J (mag) = 


H (mag)'' 


Ks (mag) = 


[3.6] (mag)'' 


[4.5] (mag)'' 


[5.8] (mag)f 


[8.0] (mag)f 


SpT'' 


Period (days)'' 


ISOY J053454.31-045413.0 


[RRS2004] NOFF W015 






14.83 


12.43 


11.29 


10.52 


10.31 


9.94 


9.26 


M4-M5 


5.1993 


ISOY J053505. 71-052354.1 


JW 276 


17.00** 


13.96=' 


12.11 


11.46 


11.16 


10.84 


10.79 






M5-M6 


20.485 


ISOY J053515. 55-052514.1 


Parenago 1872 


13.83=' 


11.55=' 


10.05 


9.11 


8.74 


8.41 


8.37 


7.50 




K4-M1'' 




ISOY J053515. 76-052309. 9 


thetal Ori E 






6.64 


6.24 


6.06 


6.94 


6.52 


6.12 




G2 IV* 


9.89520 


ISOY J053518. 03-052205. 4 


[H97b] 9209 


19.06=' 


15.30=' 




10.36 


9.30 


8.98 


8.67 






K0-K3 


5.6175 


ISOY J053526. 88-044730. 7 


Parenago 2017 


12.17= 


11.06= 


10.16 


9.59 


9.42 


9.36 


9.32 


9.30 


9.26 


K0-K2 


3.905625 


ISOY J053446. 01-044922.1' 


V1448 Ori 


17.43= 


15.61= 


13.88 


13.26 


12.95 


12.73 


12.66 


12.91 




K5= 


0.5424125 


ISOY J053605. 95-050041. 2' 


2MASS J05360595-0500413 






14.89 


14.09 


13.73 


13.51 


13.49 


13.20 




<G 


3.570535 



"Names are composed of an acronym, ISOY (Initial Spitzer Orion YSO), followed by the coordinates of the source. 

''Most common names from the literature searchable in SIMBAD database: [H97b] : iHillenbrandl 119971) . JW: I Jones fc Walkerl i 19881 ). [RRS2004] : iRamirez et al.l J2004I) . Parenago: iParenagj i 19541 ) 

= from lRebulil J2001I) 

•^from lHillenbrandl jl997l) 

"from [Carpenter et al.l J2001I) 

^frorn Megeath et al. 2012 in prep. 

gfrom lCostero et all J2006I) 

''from this work 

00 

'Probable field EBs 
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Table 2. Summary of time-series photometric observations 



Source 


Data Source^ 


Band 


Date 


range 


# epochs 


low 1 juoo'y:0'y:.oi-U'y:0'y:io.u 




[d.bj, [4.0J 


9*^ Off 9nnQ 
zo wet zuuy - 


•^n KTmr 901 1 


1 fi9 1 fi9 




Ilkirt /WFC AM 


J 


20 Ort 2009 - 


- 1 5 Dec 201 


54 








9*^ Orf 9nnQ 


Q npp 9ni n 


1 fiR 1 Rfi 

±UU, -LUU 




/^r w/ 0± 


T 




'\C\ Nm/ 901 n 
- OU IN UV iUiU 


22 




VVriy' IJJOWZ.IIII 


T 


91 Mri-^r 901 n 


9Q NTz-iir 901 


Q 




LOWELL 


^ c 


24 Oct 2009 


- 1 Dec 2009 


9 




USNO 


r 

J c 


■^0 Oct 9010 - 


- 98 Dec 9010 


2387 




SMMV 




13 Dec 1994 - 


- 26 Nov 2002 


298 




MP 


r 

J c 


1 7 Nnv 2004 


- 8 Dec 200fi 


1418 


ISOY J053515. 55-052514.1 


Spitzer/IRAC 


[3.6], [4.5] 


23 Oct 2009 


- 9 Dec 2010 


166, 166 




Ukirt/WFCAM 


J 


20 Oct 2009 ■ 


- 15 Dec 2010 


44 




APO/31" 


/c 


24 Oct 2009 - 


- 14 Dec 2010 


69 




LOWELL 


/c 


24 Oct 2009 


- 3 Dec 2009 


21 




USNO 


Ic 


30 Oct 2010 ■ 


- 28 Dec 2010 


1933 


ei Ori E 


Spitzer/IRAC 


[3.6], [4.5] 


23 Oct 2009 


- 9 Dec 2010 


66, 166 


ISOY J053518. 03-052205. 4 


Spitzer/IRAC 


[3.6], [4.5] 


23 Oct 2009 


- 9 Dec 2010 


166, 137 




WFI 


Ic 


21 Nov 2010 ■ 


- 28 Nov 2010 


7 


ISOY J053526. 88-044730.7 


Spitzer/IRAC 


[3.6], [4.5] 


23 Oct 2009 - 


- 27 Apr 2011 


108, 2821 




Ukirt/WFCAM 


J 


20 Oct 2009 ■ 


- 15 Dec 2010 


54 




USNO 




29 Oct 2010 ■ 


- 25 Feb 2011 


14927^= 




ROl 


la 


10 Dec 1995 


- 3 Feb 1997 


79 


ISOY J053446. 01-044922. 1'' 


Spitzer/IRAC 


[3.6], [4.5] 


23 Oct 2009 


- 9 Dec 2010 


92, 92 




Ukirt/WFCAM 


J 


20 Oct 2009 ■ 


- 15 Dec 2010 


54 




USNO 


/c 


29 Oct 2010 ■ 


- 25 Feb 2011 


1357 




ROl 


Ic 


10 Dec 1995 


- 3 Feb 1997 


79 




SMMV 


Ic 


13 Dec 1994 ■ 


- 22 Dec 1994 


126 


ISOY.J053605. 95-050041. 2'' 


Spitzer/IRAC 


[3.6], [4.5] 


23 Oct 2009 - 


■ 30 Nov 2011 


178, 178 




Ukirt/WFCAM 


J 


20 Oct 2009 ■ 


- 15 Dec 2010 


54 




SMMV 


Ic 


13 Dec 1994 ■ 


- 22 Dee 1994 


89 



''ROl; iRebuIll ll200lh : SMMV: [stassun et al] l|l999l '): MP: Monitor project, private communication, 
t' Probable field EEs 

■^The USNO data for ISOY J0535-0447 was binned every 9 datapoints yielding a total of 1624 datapoints 
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Table 3. Radial velocity for ISOY J0535-0447 



HJD 


RV (km/s) 


eRV (km/s) 


Instrument /Telescope 


2455463.989 


30. 


2. 


ECHL/KPN04m 


2455522.83089 


23.5 


2. 


NIRSPEC/Keck II 


2455523.06103 


23.5 


2. 


NIRSPEC/Keck II 


2455527.05545 


35.3 


2. 


NIRSPEC/Keck II 


2455543.903 


30.6 


2. 


HIRES/Keck I 


2455584.93613 


27.3 


2. 


NIRSPEC/Keck II 


2455597.9375 


38. 


2. 


NIRSPEC/Keck II 


2455635.726 


28.0 


2. 


HIRES/Keck I 



^this epoch was not used for the fit 



Table 4. Orbital parameters of 6^ Ori E 



Parameter 



Value 



Period 9.89520 ± 0.0007 

HJDo 2453285.98828'' 

asini 33.046 ± 0.106 R© 

i 73.7 ± 0.9 deg 

a 34.430 ± 0.193 Rq 

M2/M1 0.9965 ± 0.0065 

Ml 2.807 ± 0.048 Mq 

M2 2.797 ± 0.048 Mq 

T2/T1 1.12 ± 0.08 

R1+R2 12.5 ± 0.6 Rq 



Costero et al.l ([2003) 



^assumed Ti = 6000 K based on 
spectral type 
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Table 5. Orbital parameters of ISOY J0535-0447 



Parameter Value 

Period 3.905625 ±0.000030 

HJDo 2455126.26 

i 88.8 ± 0.9 deg 

a 10.0 Rq^ 

V^ 30.4 km/s 

q = M2/M1 0.06 

Ml 0.83 Mq 

M2 0.05 Mq 

T2/T1 0.55 ± 0.03 b 

Ri + R2 2.87 ± 0.03 R0 



^The semi-major axis was fixed, 
and therefore the mass ratio and the 
individual masses are estimated 

''Ti = 5150 was adopted from KO 
spectral type 



Table 6. 6^ Ori E time series at [3.6] and [4.5]. (Full table available in the electronic 

version of this article.) 



HJD 


Filter 


Mag 


Error 


2455128.92141 


IRACl 


6.793 


0.027 


2455129.30254 


IRACl 


6.802 


0.026 


2455130.43992 


IRACl 


6.807 


0.027 


2455131.24015 


IRACl 


6.828 


0.030 


2455132.78887 


IRACl 


6.803 


0.028 


2455133.3197 


IRACl 


6.806 


0.027 


2455135.00469 


IRACl 


6.769 


0.030 


2455135.67994 


IRACl 


6.779 


0.030 


2455137.19755 


IRACl 


6.776 


0.027 


2455140.00532 


IRACl 


6.785 


0.028 
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Table 7. ISOY J0535-0447 time series at [3.6], [4.5], J, and Ic bands. (Full table available 

in the electronic version of this article.) 



HJD Filter Mag Error 



2455128.00636 


IRACl 


9.273 


0.003 


2455128.26383 


IRACl 


9.309 


0.003 


2455128.50885 


IRACl 


9.261 


0.003 


2455128.87874 


IRACl 


9.267 


0.003 


2455129.25986 


IRACl 


9.288 


0.003 


2455129.75754 


IRACl 


9.204 


0.003 


2455130.39723 


IRACl 


9.245 


0.003 


2455131.19746 


IRACl 


9.258 


0.003 


2455131.85941 


IRACl 


9.258 


0.003 


2455132.08532 


IRACl 


9.300 


0.003 



Table 8. ISOY J0535-0522 time series at [3.6], [4.5], and Ic bands. (Full table available in 

the electronic version of this article.) 



HJD 


Filter 


Mag 


Error 


2455128.04907 


IRACl 


8.798 


0.092 


2455128.30655 


IRACl 


8.827 


0.093 


2455128.55158 


IRACl 


8.845 


0.097 


2455128.92141 


IRACl 


8.926 


0.101 


2455129.30254 


IRACl 


8.807 


0.091 


2455129.80022 


IRACl 


8.826 


0.092 


2455130.43992 


IRACl 


8.828 


0.097 


2455131.24015 


IRACl 


8.83 


0.097 


2455131.90210 


IRACl 


8.822 


0.099 


2455132.12801 


IRACl 


8.822 


0.096 



-48 - 



Table 9. ISOY J0535-0523 time series at [3.6], [4.5], and Ic bands. (Full table available in 

the electronic version of this article.) 



HJD 


Filter 


Mag 


Error 


2455128.04658 


IRACl 


10.672 


0.044 


2455128.30406 


IRACl 


10.671 


0.044 


2455128.54909 


IRACl 


10.653 


0.042 


2455128.91892 


IRACl 


10.681 


0.043 


2455129.30005 


IRACl 


10.669 


0.044 


2455129.79774 


IRACl 


10.650 


0.045 


2455130.43743 


IRACl 


10.688 


0.045 


2455131.23766 


IRACl 


10.662 


0.045 


2455131.89961 


IRACl 


10.667 


0.044 


2455132.12552 


IRACl 


10.658 


0.042 



Table 10. ISOY J0534-0454 time series at [3.6], [4.5], and Ic bands. (Full table available 

in the electronic version of this article.) 



HJD 


Filter 


Mag 


Error 


2455128.00073 


IRACl 


10.472 


0.039 


2455128.25819 


IRACl 


10.462 


0.038 


2455128.50321 


IRACl 


10.471 


0.039 


2455128.87311 


IRACl 


10.451 


0.037 


2455129.25424 


IRACl 


10.469 


0.039 


2455129.75191 


IRACl 


10.452 


0.039 


2455130.39161 


IRACl 


10.473 


0.036 


2455131.19183 


IRACl 


10.480 


0.039 


2455131.85379 


IRACl 


10.466 


0.039 


2455132.07970 


IRACl 


10.460 


0.037 
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Table 11. ISOY J0535-0525 time series at [3.6], [4.5], and Ic bands. (Full table available 

in the electronic version of this article.) 



HJD Filter Mag Error 



2455128.04668 


IRACl 


8.337 


0.033 


2455128.30416 


IRACl 


8.358 


0.032 


2455128.54919 


IRACl 


8.380 


0.033 


2455128.91902 


IRACl 


8.432 


0.031 


2455129.30015 


IRACl 


8.456 


0.037 


2455129.79818 


IRACl 


8.494 


0.036 


2455130.43769 


IRACl 


8.498 


0.036 


2455131.23792 


IRACl 


8.443 


0.037 


2455131.89987 


IRACl 


8.389 


0.035 


2455132.12578 


IRACl 


8.365 


0.033 



